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ABSTRACT 

This  report  is  concerned  with  the  theoretical  evaluation  of  the 
buckling  strength  of  a  truncated  hemisphere  under  axial  tensile  load. 
The  edges  of  the  shell  are  assumed  to  be  restrained  from  moving 
radially  or  from  rotating.  Theoretical  results  were  obtained  by 
Vlasov'  s  small  deflection  theory  and  Galer  kin's  method.  The  quick 
convergence  of  the  series  solution  is  demonstrated.  Comparison  of 
theoretical  values  with  a  few  available  experimental  results  is  given. 
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SYMBOLS 

coefficients  of  series  expansion  for  radial  deflection  and 
stress  function,  respectively 

mean  radius  (see  Figure  2) 

Young's  modulus 

sin  a 

number  of  buckle  waves  in  the  circumferential,  direction 

axial  tensile  load  per  unit  circumferential  length  at  the 
shell  equator 

radial  load  per  unit  middle  surface  area 
additional  force  components  in  the  buckled  shell 
membrane  force  components  prior  to  buckling 
shell  thickness 

radial  buckling  displacement  of  a  middle  surface  point 
sin  $ 

altitude  angle  of  the  truncated  edge  (see  Figure  2) 
spherical  coordinates  (see  Figure  2) 

Poisson's  ratio 
stress  function 
Laplace's  operator 

(ja’2/a 2ftp2)  -tan  p  (3/a 23p)  +  sec2  <f>  (32/a2392)J 
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I.  INTRODUCTION 

A  truncated  hemisphere  may  buckle  under  an  axial  tensile  load  that  exceeds 
a  certain  critical  value  (Figure  1)  because  of  compressive  hoop  stresses.  The 
buckled  shape  is  similar  to  that  of  a  cylinder  buckled  under  radial  pressure, 
with  one  half  wave  in  the  axial  direction  and  many  small  waves  in  the  circum¬ 
ferential  direction.  The  present  theoretical  study  treats  the  case  where  both 
the  truncated  edge  and  the  edge  at  the  equator  are  restrained  from  moving 
radially  or  from  rotating. 

In  the  analytic  work,  Vlasov's  equations  of  equilibrium  and  compatibility 
condition  (Reference  1)  are  used.  These  equations  are  satisfied  by  the  use  of 
Galerkin's  method  (Reference  2)  in  conjunction  with  expressions  for  the  radial 
displacement  and  the  stress  functions  which  satisfy  the  approximate  boundary 
conditions . 

Numerical  results  of  theoretical  buckling  loads  for  shells  of  various  geo¬ 
metric  configurations  are  given  by  curves  as  well  as  by  simplified  algebraic 
formulas,  and  compared  with  a  few  experimental  data. 

II.  THEORY 


Using  spherical  coordinate  system  0  and  8,  as  shown  in  Figure  2,  the 
equilibrium  condition  in  the  radial  direction  for  a  spherical  shell  element  can 
be  expressed  by  the  following  equation: 


where 


1 

12  (1-p2) 


4_4 
a  V 


2  2 
a  V 


a2  ^  a  2  a2 

— t  -  tan  0  i—  +  sec  0  — 

a/  **  ae 


7  = 


2  2 
7  7 


(1) 


and  $  is  the  stress  function,  in  terms  of  which  the  additional  force  components 
in  the  buckled  shell  wall  are  given  by 
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Fig.  1.  Truncated  Hemisphere  Buckled  by  Axial  Tension 


Equations  (2)  approximately  satisfy  the  equations  of  equilibrium  for  a  shell 
element  in  both  the  meridional  and  the  circumferential  directions.  The  stress 
function  ■§  arid  the  radial  deflection  w  are  also  related  by  the  compatibility 
equation 


Equations  (1)  and  (3)  are  Vlasov's  equations  (Reference  1,),*  expressed  in  spherical 
coordinates . 

In  dealing  with  the  problem  of  stability  of  a  spherical  shell,  we  must  take 
into  account  the  radial  component  of  stress  existing  prior  to  buckling.  With  the 
assumption  that  these  prebuckling  Stresses  are  adequately  represented  by  the 
membrane  state  of  stress,  we  have 


T*o  =  -  TS0  =  P  *  <4> 

During  buckling,  because  of  curvature  changes,  these  finite  membrane  forces 
Contribute  a  radial  component  equal  to 

(5) 


q  = 


-  sec 


.2 

*  §  w  x 
0  — T  + 

a© 


tan 


sec  <p 


which  is  used,  in  Equation  (T). 


Page  5 


The  boundary  conditions  chosen  to  be  satisfied  by  the  radial  deflection  w 
and  the  stress  function  $  are  those  corresponding  to  shell  edges  that  are 
restrained  from  moving  radially  or  from  rotating,  motion  in  the  meridional 
and  circumferential  directions  being  restrained  insofar  as  rigid  body  motion  is 
concerned.  These  conditions  may  be  expressed  as 

w  =  $  =  Tp  =  v  =  0  at  * =  0,Q  •  m 

or,  for  Equation  (2),  as 

w  =  -  j  =  ^  =0  at  p  ,=  0,  a  .  (6b) 


A  possible  method  of  solving  the  problem  would  be  to  first  assume  an 
arbitrary  radial  deflection  function  that  satisfies  the  radial  deflection  boundary 
conditions  (6);  then  Equation  (3)  could  be  solved  for  the  stress  function  f  in 
terms  of  a  particular  solution  involving  w  and  the  general  solution  of  the  homo¬ 
geneous  equation,  the  constants  of  integration  being  determined  by  the  stress 
function  boundary  conditions .  Finally,  the  radial  deflection  function  and  the 
derived  stress  function  can  be  substituted  into  Equation  (I),  which  then  could  be 
solved  by  means  of  the  Galerkin  method  (Reference  2).  Because  the  solution  for 
the  stress  function  in  terms  of  the  radial  deflection  function  is  difficult,  however, 
this  method  of  solution  can  be  replaced  by  the  equivalent  process  of  choosing 
arbitrary  functions  which  satisfy  the  appropriate  boundary  conditions  for  both 
the  stress  function  t  and  the  radial  deflection  function  w  and  solving  both 
Equations  (1)  and  (3)  by  the  Galerkin  method.  Thus, 


ify*"***- 


a  cos  $  d0  dp  =  0 


(7a) 
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}2tt 


2 

a  cos  0  46  d^  ••,=  0 


where  w  and  $  are  subjected  to  the  boundary  conditions-  given  by 
For  convenience,  the  following,  substitutions  are  introduced: 


X  = 

sin  <p  ,  l 

=  sin  a  , 

.a 

1  2  ; 

b<p 

T1  -  x  * 

|  ,  etc.  , 

2-2 

2.  a2 

„  a 

a  V  =  ! 

(i-x  )-% 

"  2x  W  *  ' 

1.  | 
i-xz 


The  boundary  conditions  (6b)  then  become 

3w 


w  = 


=  0  at  x  =  0,  i 


at 


5  =  lx  =  0  at  x  =  °«  1 

Suitable  expressions  for  w  and  i  that  satisfy  conditions  (9)  may  be  taken  as 


w 

T 


N 


A  x11  (x-jt)n  cos  m9 
n=2,  3,  .  .  .  n. 


t 


B  xn  (x-X)n  cos  m-0 
n 


Ea  t  ii=2,  3,  .  .  . 

The  substitution  of  Equations  ( 10)  into  Equation  (7)  then  leads  to 


f!  (A  e  .  i  B  d  .1  =  0 
^  l  “  nj  n  ny) 


n=2 


/l“ 

'f  d  .-If ' 

A  -  ; 

B  e  . 

lll2 

u-nM 

aj  nj  Et  nj 

n 

n  nj  j 

=  0 


(ID 


j  =  2,  3,  4,  .  .  .N 


^fejjr2)  i2n+2j'3  {^T  (n+j )'.  (ri+j-4)!  -  4  (n+j*.  1)1  (n+j -3)1 


+  3^  (n+j- 1)1  (n+j-1)'] 


'  (2n^2j'-V).  i2n+2j'1  |(7n2-5n-2)(n+j)I  (n+j -2)! 

z 

-  kn  (n+1 )  +  (6+2m2)j  |tt+j-l)lj  -  4  (n2-n-2)(n+j+l)l  (n+j -3)'. 


+  (n- 1  )(n  -  2  )(n  -  3 )( n+j +2 )  1  (n+j -4)1 


+  WtFTTT  *2ti+2j+1  +  W 


-  (4  n3+12n2-14n+4)(n+j+l)l  (rt+j-1)!  +  (6n2+18n+l4)(n-l)(n+j+2)! 


-  4  (n- 1  )(n  -4)(n+ j+3 )  1  (n+j-3)'  +  (n-l)(n-2)(n-3)(n+j+4)l  (n+j-4)!j 

.  /o__4  ^ _ 2^  (n+j)!  (n+j+2)!  <2n+2j+3  A/~  4  , «,  -2..  (n+j+4)l,(n+j);  ... 

+  (2m  -6m  )  ••  (fc+ffi?)!--  1  H3m  ' 10rn  ’  (2n+2j+5)l~^ 1 


2n 

(2ti+2j 


_je2n+2j-l|(n_i)(n+j)l  (n+j.2),  -  n  §ii+jTl}lj  } 
_fc2n+2j+l  |n+i+  n^  |n+j)!]  2  .  2(n+l)(n+j+l)l 


+  (n-l)(n+j+2)l  (n+j 


:j-2)'-| 


_ 2  (n+j+2)'.  (n+j)'.  ,2n+2j+3  _2  (n+j)*  (n+j+4)'.  #2n+2j+5 

"  (2n+2j+3)t  1  m  (2n+2j+&)!  * 


For  a  nontrivial  solution  for  A  and  3B  ,  the  determinant  of  their  coefficients 

n  n 

in  Equations  (11)  must  vanish,  yielding  the  stability  criterion  for  the  problem. 
The. lowest  root  of  the  determinal  equation  is  the  critical  value  of  the  load 


III.  NUMERICAL  RESULTS  AND  DISCUSSION 

Computations  for  the  buckling  tension. load  were  carried  out  for  shells  with 
a  .=  20,  30  and  40  degrees  and  a/t  ranging  from  200  to  1600.  Four  terms  each 
were  used  in  the  summation  of  Equations  (10).  The  output  is  given  by  the  non- 
dimensional  quantity  P/Et  cos2  a  against  a/t  as  shown  in  Figure  3.  It  may  be 
seen  that  within  the  range  of  a  (that  is,  a  is  less  than  40  degrees),  the  buckling 
strength  of.  the  shell  is  far  less  influenced  by  a  than  by  a/t.  The  curve  shown 
in  Figure  3  may  also  be  expressed  by  the  following  approsdmate  formula: 

2  m1-12 

P/Et  Cos6  o  .=  1.57  fg  (15) 


In  the  limiting  case  where  the  radius  of  the  sphere  approaches  infinity,  the 
problem  is  reduced  to  a  narrow,  infinitely  long  strip  which  is.  subjected  to  com¬ 
pression  along  the  length  and  tension  across  the  width  of  the  strip.  In  this  case, 
the  buckling  stress  is  of  the  form  (Reference  3,  p.  337) 


where  k  is  a  constant  coefficient. 


100  200  300  400  600  800  1000  2000 

q/t 


Fig.  3.  Buckling  Stress  of  Truncated  Hemisphere  (p  =  0.  3). 


We  recall  (Reference  3,  pp.  450,  496)  that  the  classical  buckling  stress  for 
Cylindrical  shells  under  radial  pressure  is 


„  .1 

4(1-  A 


and  that  the  classical  buckling  stress  for  spherical  shells  under  external  pres¬ 
sure  is 


We  may  see  from  Equations  (4)  and  (15)  that  the  buckling  characteristic  of  the 
present  problem  is  similar  to  the  sphere  problem  rather  than  the  cylinder 


Test  data  for  five  samples  of  various  geometry  were  furnished  by  Douglas 
Aircraft  Company  and  are  reproduced  here  in  Table  1.  The  theoretical  buckling 
load  parameter  P/ Et  obtained  by  using  one  term  in  the  summation  of  Equation 
(10)  as. the  first  approximation  and  succeedingly  using  more  terms  in  the  sum¬ 
mation  as  higher  appr oximation s  is  also  shown.  We  note  that  the  convergence 
of.  the  series  solution  is  good  and  that  using  four  terms  in  the  summation  is 
adequate  for  practical  purposes.  The  theoretical  values  are  found  to  be  about 
twice  greater  than  the  experimental  value s ,  a  difference  rather  common  in 
buckling  problems.  In  Table  1,  also  given  are  the  theoretical  cir cumf e r ential 
buckling  wave  numbers  associated  with  the  found  theoretical  minimum  buckling 
loads.  These  numbers  seem  quite  high  in  comparison  with  the  wave  numbers 
found  by  test.  However,  a  study  of  the  curves  for  the  theoretical  buckling  load 
versus  wave  numbers  in  Figure  4  leads  to  the  belief  that  the  actual  buckling 
wave  number  probably  is  determined  by;  the  initial  imperfection  of  the  shell 
since  the  curve  in  the  neighborhood  of  the  theoretical  minimum  point  is  rather 
flat. 

The  subject  problem  is  an  interesting  one  as  it  is  a  particular  example 
that  an  elastic  system  is  buckled  by  tension,  although  the  real  cause  of  the 
instability  is  the  compressive  hoop  stress.  Certainly,  there  remain  many 


areas  for  future  studies,  among  them  determination  of  the  buckling  load  by  large 
deflection  theory  and  consideration  of  initial  imperfections.  The  testing  results 
presented  here  are  rather  scarce  and  scattered,  hardly  enough  to  give  any  con¬ 
clusive  indication.  Extensive  and  systematic  experiments  in  this  respect  are 
recommended. 
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